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ABSTRACT 


Computer simulation studies of a frequency domain adaptive beamforming 
algorithm are presented. These simulation studies were conducted to determine 
the multiple broadband target localization capability and the full angular coverage 
capability of the algorithm. The algorithm was evaluated at several signal-to-noise 
ratios with varying sampling rates. The number of iterations that the adaptive 
algorithm took to reach a minimum estimation error was determined. Results of 
the simulation studies indicate that the algorithm can localize multiple broadband 


targets and has full angular coverage capability. 
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I. INTRODUCTION 


This thesis is but a part of an ongoing research project to develop new sonar 
signal processing algorithms capable of quickly and accurately solving target lo- 
calization problems. Present technology and doctrine dictate that several lines of 
bearing to a target be obtained before a sonar fire control solution can be com- 
puted. Obtaining these lines of bearing is a time consuming and often dangerous 
task due to the increased probability of counterdetection and, as a result, evasive 
maneuvering and defensive action on the part of the target. A sonar system ca- 
pable of providing timely, accurate target localization while minimizing own ship 
maneuvering would result in longer firing ranges and, therefore, a reduction of the 
threat to one’s own ship. 

Several recent papers [Refs. 1-3] have discussed the application of a complex 
least-mean-square (LMS) adaptive algorithm [Ref. 4] to bearing estimation prob- 
lems using a linear array of sensors. When a linear array is used, only a bearing 
angle ~, to the source can be estimated, an estimate of the depression angle 6, can- 
not (Figure 1.1). While an estimate of the bearing angle is useful in the localization 
problem, a better tool would be an algorithm which provides estimates of both the 
bearing angle and the depression angle. A frequency domain adaptive beamform- 
ing algorithm for planar arrays that solves multiple, broadband target localization 
problems and provides estimates of both the bearing angles and the depression 
angles from the center of the planar array to the targets has been developed and 


initially tested by Ziomek and Chan [Ref. 5]. 





Figure 1.1. Target Location Geometry. 


The research performed in this thesis was to continue the work of Ziomek 
and Chan [Refs. 5, 6] and fully evaluate the capabilities of the complex, least- 
mean-square, frequency domain, adaptive beamforming algorithm they developed 
via computer simulation studies. The computer simulation studies were designed 
to test the algorithm’s multiple broadband target localization capability, its full 
angular coverage capability, and its angular resolution as a function of the input 
signal-to-noise ratio (SNR) at a single element in the array, sampling rate, harmonic 
number, and the number of iterations of the algorithm. 

Each target was modeled as a broadband sound source. As a result, the 
frequency spectrum of the output signal from each element of the planar array 
contains several frequency components. An estimate of the bearing and depression 
angle for each frequency component is provided as a result of processing the output 
frequency spectrum from each element in the array through the frequency domain 
adaptive beamforming algorithm. Therefore, if each target exhibits at least one 
unique frequency component (or spectral line), then all targets can be located. 

Full angular coverage is the ability to localize a target regardless of its relative 
position to the array. The broadside case is the easiest since it is at this position 
that the far-field beam pattern beamwidth 1s its narrowest. The endfire case is the 
most difficult since the far-field beam pattern beamwidth is the broadest at this 
point. The full angular coverage and multiple broadband target capabilities were 
tested simultaneously. 

Baseline results were the first assembled. Baseline results are defined as bear- 
ing and depression angle estimation errors (measured in degrees) as a function of 
sampling rate, harmonic number, and the number of iterations of the algorithm 


for the “no noise” case. Identical cases were run using additive, zero mean, white, 


guassian noise to corrupt the output signals from each element of the planar ar- 
ray. Average bearing and depression angle estimation errors were then plotted as 
a function of the input SNR at a single element of the planar array, the sampling 
rate, and the harmonic number. 

Chapter II describes the theory used in the development of the frequency do- 
main adaptive beamforming algorithm. The construction of the direction cosine 
estimates, the angle estimates, and phase “unwrapping”, integral parts of the al- 
gorithm, will be presented in detail. 

Chapter III contains computer simulation results and an explanation of these 
results. The results which are presented consist of graphical representations of 
the average estimation errors of the bearing and depression angles at two distinct 
signal-to-noise ratio levels for four cases. These cases include targets evaluated 
at broadside, endfire, random placement, and targets which share a spectral line. 
Conclusions concerning the effect of harmonic number, sampling rate, number of 
iterations of the algorithm, and SNR are made. The Appendix contains tabular 
numerical data for all results. 

Chapter IV will present final conclusions and recommendations for further 
research. 

Results of this thesis research work have been incorporated into a research 
paper by Ziomek and Behrle which was presented at the Twenty-First Annual 


Asilomar Conference on Signals, Systems, and Computers [Ref. 7]. 


II. THEORY 

This chapter is designed to present a brief, yet comprehensive, development of 
the frequency domain adaptive beamforming algorithm which is evaluated in this 
thesis. The goal of the algorithm is to process the output electrical signals from 
a planar array of sensors and provide estimates of direction (both bearing and 
depression angles) and frequency content of the acoustic fields incident upon the 
array. The development of the algorithm presented here is similar to the analysis 
section of a paper by Ziomek and Behrle (Ref. 7], but this chapter provides a more 
detailed explanation for several parts of the algorithm. A more in-depth analysis 
of the algorithm is provided by Ziomek and Chan [Ref. 5] and Chan [Ref. 6]. 

The frequency domain adaptive beamforming algorithm is based on using the 
complex frequency domain data R(q,m,n) from all M x N elements of a planar 
array. The objective of the adaptive filter used in the algorithm is for the filter 
to converge to a set of phase weights such that the array output signal will match 
a reference signal. From this set of phase weights, estimates of the bearing and 
depression angles as a function of the harmonic number g can be made. 

Consider a M x N planar array of point source elements lying in the XY 
plane (Figure 2.1). These elements are equally spaced where d, and d, are the 
interelement spacings in the X and Y directions, respectively, and M and WN are 
the total number (odd) of elements in the X and Y directions, respectively. The 
random baseband output electrical signal, r(@, m,n), is composed of a deterministic 
signal, y(@,m,n), and a random receiver noise component, n(£,m,n). The output 


signal at time instant 2 and element (m,n) in the array is given by 





Figure 2.1. Planar Array Geometry. 


where 
J 2 oe (222) 
M' =(M —-1)/2, (2.3) 
N' =(N —1)/2, (2.4) 
and 
L>2K+1 (2:5) 


is the total number of time samples that must be taken per element to avoid aliasing 
when the deterministic signal y(@,m,n) is composed of K harmonics [Ref. 8: p. 
164]. The minimum value for L (i.e, L = 2K +1) is obtained by sampling at 
the Nyquist Rate. If T, seconds is the length of the data record (or fundamental 
period) recorded at each element of the array, then the fundamental frequency (or 


the FFT bin spacing) is given by 
fo =1/To Hz. (2.6) 


Therefore, the highest frequency which can be contained within the deterministic 


part of the received signal is 

fmaz = Kfo Hz (2.0) 
and the minimum sampling frequency of the received output electrical signal r(2, m,n) 
is 


fs = L/Ty Samples/sec (2.8) 


where L must satisfy equation 2.5 [Ref. 8: p. 164]. 

To obtain the complex frequency domain samples of the received signal, the 
discrete Fourier transform (DFT) with respect to the time index @ of equation (2.1) 
is taken. This action yields 

R(q,m,n) = Ylq,min) + NG.) = ee ees 
m= —M",...,0,...,M’ (2.9) 
n=-—N',...,0,...,N’ 
where the index g represents the harmonic number. 


rng 


Consider a single general plane-wave field, g (: i - ), propagating in the 





+n direction, incident upon the planar array as shown in Figure 2.2. If ny isa 
unit vector, g(t) an arbitrary baseband function, and c the speed of sound in the 
medium measured in meters per second, then the deterministic part of the output 
electrical signal at time instant @ and element (m,n) in the array is given by [Ref. 


8: p. 160} 


ugmdx + vondy (2 10) 
a : 


y(@,m,n) = (rs + 
where T's is the sampling period in seconds, and wp and vo are dimensionless 
direction cosines with respect to the X and Y axes, respectively. 

Upon taking the DFT of equation 2.10, we obtain the corresponding frequency 


spectrum given by [Ref. 8: pp. 162-166] 
NAC hag) — le, exp(+ j27q fouomdx /c) exp(+ j27q fovondy /c) (Zoi 


where cy, g = —L’,...,0,...,£" are the complex Fourier series coefficients that 
can be used to represent the baseband function g(t) by a finite Fourier series with 


Kk harmonics during the time interval |t| > Ty) /2. The coefficients are given by 


€qi = 1090.0) eee ee (2.17) 


>) 
» 


Figure 2.2. General Plane-Wave Field g (+3 ri ) Propagating in the 
tno direction. 


The dimensionless direction cosines, ug and vg, are given by 


ug = sin A COS Wo (2a 


Vp = sin Bp sin yo (2.14) 


where @ is the depression angle and wo is the bearing angle. These angles are not 
known a priori. 

If we consider multiple plane waves incident upon the planar array, equation 
(2.10) can be generalized and the deterministic part of the output electrical signal 


at time instant @ and element (m,n) is 


(2.15) 


uormdx + ante | 
Cc 


WG Oy 0) S- gk (ers i 
k 
where gx(t), vox; and vox are the arbitrary baseband function and direction cosines 
associated with the kth sound source. 
A brief description of the development of the frequency domain adaptive beam- 
forming algorithm is presented next. The discussion that follows is based, in part, 
on the material in Ref. 7. 


Define the complex estimation error at harmonic gq as 
A 2 
e(q) = s(q) — S(q). (2.16) 


The reference signal s(q) is defined as 


M' N‘ 
(=~ YY [Ramn)lexp(+i2R(a,0,0)) (2.17) 
m=-M'n=—-N' 
and 
1 M' N’ 
()=TW dD dL am)d(g,n)R(q,m,n) (2.18) 
m=—M'n=-N' 
= c’ (q)R(q)d(q)/(LMN) 
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is the estimate of s(q), where c(q,m) and d(q,n) are the unit magnitude complex 


weights in the .Y and Y directions, respectively, c(q) and d(q) are the M x 1 and 


N x 1 complex weight vectors in the X and Y directions, respectively, which are 


given by 


c(q) = [e(q, -M"),...,¢(q,0),-.-,e(g, M’)]* 


and 
d(q) = [d(q,-N’),...,d(q,0),...,d(q, N’)]” 


and R(q) is the M x N complex data matrix given by 


R(q,-M',-N’) ... R(q,-M',0) ... R(q,-M',N') 
Pag) = R(q,0,-—N') eee R(q, 0,0) ae R(q,0, N‘) 
R(q,M',-N')  ... R(q,M',0) ... R(q,M’',N’') 


Now define the (M + N) x 1 complex weight vector w(q) as follows: 


0? Li] 


We can obtain the original weight vector in the X direction via 


c(q) = A w(q) 


where 
I | 0 


She SM MxwN 


(2.21) 


CRey) 


isa M x(M +N) matrix, J isa M x M identity matrix, and 0 isa M x N null 


matrix. Similarly, we can obtain the original weight vector in the Y direction via 


d(q) = B w(q) 


where 


ILL 


(2.25) 


0 aia i 


DNA es (2.26) 


is a N x (M+ N) matrix, 0 is a N x M null matrix, and J is a N x N identity 
matrix as indicated. 


Substituting equations (2.23) and (2.25) into equation (2.18) yields 


$(q) = w’ (q)Z(q)w(9)/(LMN) (2.27) 


where 


Z(q) =A’ R(q)B (2.28) 


isa (M+WN)x(M+WN) complex matrix. The complex weight vector that minimizes 


the mean-square error E-{|e(q)|*} is given by [Ref. 5, Ref. 7] 


Wi4) (q) a w;(q) AR 2;e;(q) |Z(q) 7 a Z7 (q)| " w; (q),% = 0,1,2,... (2.29) 


where 


e:(q) = s(q) — §:(g) (2.30) 


$;(q) = wi (q)Z(q)w,(¢)/(LMN) (2.31) 


is the estimate of s(q) after the it* iteration, and the step-size parameter p; is 


given by 


TES Vy (05 eos) 4 SSO ee. (2732) 


where po is a constant, and is equal to the inverse of the sum of the signal and 
noise power, af and a7, respectively, at the center element of the array. After each 
iteration of the algorithm, each component of the complex weight vector w,.,(q) 


is normalized by its respective magnitude in order to maintain unit magnitude. 
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Once the complex weight vector, w,,,(q), converges to a steady-state value, 
w,,(q), the steady-state complex weight vectors c,,(q) and d,,(q) can be obtained 
from equations (2.23) and (2.25), respectively. The estimates of the depression and 
bearing angles at harmonic gq, 9o(q) and wo(q), respectively, are given by [Ref. 8: 
p. 175] 


1/2 


69(q) = sin” { (ah S(q)]° + [es5(a)]" } | gel (2.33) 
and 


vo(q) = tan [6f°(¢)/ag>(q)], ¢ #0 (2.34) 


where 7*°(q) and 64°(q) represent estimates of the direction cosines obtained by 
using a least-squares fit to the “unwrapped” steady-state phase weights 6%,(q,m) 


and 0% .(q,n), respectively. In the absence of noise 


u —L’,. yulee 
6. .(q,m) = t2rqfouo(q)mdx/c, M’.....0....M (2935) 
and 
u =|) Deer 0 een Oy 
$39(4,2) = £2rqfovo(q)ndy/c, ee 0 IN (2.36) 


It should be noted that if r(@,m,n), y(@,m,n), and n(£,m,n) of equation 
(2.1) are baseband complex envelopes, then qfp must be replaced by (f, + qfo) in 
equations (2.35) and (2.36), where f, is the carrier frequency in hertz. In addition, 
equations (2.33) and (2.34) must now be evaluated at gq = 0. Substitution of 
equations (2.17) and (2.18) into equation (2.16) yields the following expression for 
the steady-state estimation error 


M' N! 
1 
el) = rape DD Ilan) expt /R(q,0,0)] 


— exp {+j{0%,(q,m) + 6%(q,n) + ZR(q,m,n)]} (2.37) 
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where 


R(q,m,n) = |R(q,m,n)|exp[+j 2 R(q,m,n)], (2.38) 
Cas(q,™m) = a55(g,m) exp[+jO3,(q,m)] = exp [+j93,(q,m)], (2.39) 

and 
dys(q,n) = bs5(q,) exp[+7 $5,(4,”)] = exp[+743,(4,7)] (2.40) 


where a,,(q,™m) = 1 and b,,(q,n) = 1 are real, unit magnitude, amplitude weights 
and 6“%(q,m) and ¢%(q,n) are real, “wrapped”, phase weights. 

To ensure obtaining the correct depression and bearing angle estimates, 9o(q) 
and bo(q), respectively, the steady-state phase weights need to be “unwrapped.” 
The “unwrapped” steady-state phase weights can take on values outside the closed 
interval [—z, 7] and therefore ensures full angular coverage capability (i.e., 0 < 
6o(q) < 7/2 and 0 < w%o(q) < 27). The use of “wrapped” phase weights is necessi- 
tated by computer programming limitations. 

The question now arises how to obtain the “unwrapped” steady-state phase 
weights in order to obtain least-squares estimates of the direction cosines at each 
harmonic. The first step is to force the “wrapped” steady-state phase weights, 
6’ (q,m) and ¢%(q,n), to be equal to zero at the center of the planar array (i.e., 
m = 0,n = 0). This is accomplished by multiplying each component of c,,(q) 


’ 


and d,,(q) (equations (2.19) and (2.20)) by exp[—76%(q,0)] and exp[—j7¢%(q, 0)], 
respectively. 


The next step is to obtain rough estimates of the direction cosines, ug(q) and 


Vo(q). The “wrapped” steady-state phase weight, 6%(q,m), can be expressed as 
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(30,5... 
Bergen cn 2 Outi ama.) clacim, q=—=—L ,...,0,...,L' (2.41) 


i= Ve ee 


where 7 is chosen to ensure that the value of 6%(q,m) is within the interval [—7, zx]. 


In the computer simulations presented in this thesis, the interelement spacing is 
d, = dy = Amin/2 (2.42) 
where the minimum wavelength is 
ee Cana (2.43) 


and 


(ela = Bf He. (2.44) 


Substituting equations (2.43) and (2.44) into equation (2.42) yields 
de diese) (2h) (2.45) 


where L’ is defined by equation (2.2) and f, is defined by equation (2.6). Substi- 


tuting equation (2.45) into equation (2.41) yields 


(at) > ol a ee 


Part Lmiuniaim 4. 2m im—_—if,...,0,...,M (2.46) 


1d 


By letting 2 = 0 in equation (2.46), the following rough estimate of direction cosine 


u, at harmonic g and element m is obtained: 
P Lge 
to(q,m) = aaa 3(9,m), g#0, m=1. (2.47) 


For purposes of this thesis, the rough estimates of u,(q) used were obtained by 
averaging the results obtained by evaluating equation (2.47) at elements m = +1, 


that is, 


tio(q) = 0.5[d,(q,1) + %.(¢,-1)], 9 #0. (2.48) 


The most difficult case for the algorithm is to locate the highest harmonic, 
that is, g = +L’, when it exists at endfire (6, = 90°) relative to the planar array. 
Let us assume that the value of the direction cosine at the highest harmonic is 


equal to one. Therefore, if g = +L’, equation (2.46) reduces to 


w . 1 =O; eee 
gv (L',m) = mz + 2ir, Vo. Osis (2.49) 


Evaluating this relationship at element m = +1 to obtain the rough estimate of 


uo(q) yields 
Oe a) See ee 0 ool bya) eee (2.50) 


and since the “wrapped” steady-state value must be in the interval [—7, z], 2 must 


be equal to zero. Equation (2.50) is reduced to 
PUG sei) = seg. (2tae 


If the plus sign is chosen in equation (2.47), and m = +1 and gq = L’, then 


evaluation of equation (2.47) yields 


io(L', +1) = 6% (L',+1)/ +7). (2.52) 
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Upon substituting equation (2.51) into equation (2.52) we obtain 
tree, 1) = | (2.53) 


which is the expected result in the absence of noise since the assumption was made 
that the value of the direction cosine was equal to one at the highest harmonic. 
Evaluating equation (2.49) at a different element (i.e., m 4 +1) would yield 2 # 0 
and, as a result, an incorrect rough estimate of uo(L’) = 1. This validates the 
use of elements m = +1 to obtain the estimate of the direction cosines (equation 
(2.47)). 

The third step is to generate rough estimates of the “unwrapped” steady- 
state phase weights, 6" (q,m), by replacing u,(q) in equation (2.46) with u,(q) 


from equation (2.48), that is 


NE ne Oa le 


ZecliGh m) a +(q/L')mi(q)x, aa =f" ie 0, ex M! cs (2.54) 
Therefore, if 
ie 3) 
in <|6°.(g,m)|< (+2), g=—L’,...,0,...,£' (2.55) 








FOS ig en lO el 
then 


a 1d) Oe 


6°.(q,m) = 0),(¢,m)+ sgn (6,(4,m)| C= hic ——. 0, ..., 0 (2.56) 


Vi 


where sgn|[ | is the sign function. Note that the right-hand side of equation (2.56) 
is simply the noise corrupted version of the right-hand side of equation (2.35). 
The method of least-squares is then used to fit a straight line to the unwrapped 
phase weights 6%,(q,m) as a function of element number m for each harmonic g. The 
phase weights 6%,(q,m) were computed using equation (2.56). The least-squares 


slope at harmonic q is given by 
M’ M' 
Srs(q)= >) mO3.(q,m)/ > m’. (2.57) 
m=-—M; m=—M! 


Substituting equation (2.45) into equation (2.35) yields 


mie eames 


Bsa(qim) = £(a/L')ruo(aym, hag ge (288) 
in the absence of noise. Therefore, since in the presence of noise 
u ms 1. q=—L', gO}. ay bf 
O53(9,™) © Srs(q)m + brs(q), _M'....0.....M (2.59) 
where 
: a we 
bs(g)= 77 dD, %Ss(a,™) (2.60) 


m=—M’' 


is the least-squares “y intercept” at harmonic g, comparing equations (2.58) and 
(2.59) yields 


Ein (Gi camera) ae 0 (2.61) 


which is the least-squares estimate of direction cosine u, at harmonic q. 
In a similar derivation, it can be shown that the least-squares estimate of the 


direction cosine v,(q) at harmonic q is given by 


+67°(q) = [L'/(q7)]Sxs(q), @ #0 (2.62) 
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where now 








NS N' 
Sus(q)= > ndi,(q,n)/ do ne? (2.63) 
n=—N! n=—wN’ 
and, if 
IO... . 
ae b% (q,n) SCS? rae 0) eee’ (2.64) 
peo... (eG 
then 


ee 
b44(q.n) = (qn) + sgn [d%(a,n)] E+ 1)t, g = —L',...,0,...,L! (2.65) 


DS IN eee a he 


where 
~ as me (1 — 2 bre | eer By 
O33(9; n) ate +(q/L )\nvo(q)r, i EINE mn 0, me _N' (2.66) 
Bo(q) = 9.5 [d0(g,1) + 4.(¢,-1)], 9g #0 (2.67) 
and 
i 
UG) — Og, (Gs 0, = 1. (2.68) 
qnT 


Using the least-squares estimates of the direction cosines #/°(q) and i/?(q), 
that is, equations (2.61) and (2.62), respectively, yield estimates of the depression 


and bearing angles, 6,(q) and ,(q), respectively, using equations (2.33) and (2.34). 
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Ill. RESULTS 


Significant computer simulation results from four unique test cases are pre- 
sented in this chapter. The test cases were designed to test the algorithm’s multiple 
broadband target localization capability, its full angular coverage capability, and 
its angular resolution as a function of the input SNR at a single element in the 
array, sampling rate, harmonic number, and the number of iterations. 


The test cases were comprised of the following: 


(1) Case 1 - a single broadband target located at broadside relative to the planar 
array; 


(2) Case 2 - three broadband targets located at random positions; 


(3) Case 3 - a single broadband target located at endfire relative to the planar 
array; and 


(4) Case 4 - three broadband targets, two of which share a common harmonic. 


The simulation results of these four test cases are based upon processing the 
output electrical signals from a 7 x 7 planar array of equally spaced hydrophones. 
The acoustic field incident upon the planar array was, in general, the sum of several 
plane-wave fields travelling in different directions. Each plane-wave field consisted 
of an arbitrary number of harmonics (spectral lines) emanating from one of the 
broadband sound sources. As a result, the output electrical signal at each element 
of the planar array was composed of an arbitrary number K of harmonics, all with 
identical amplitudes of unity. 

The fundamental frequency for all test cases was chosen to be f, = 1000 
Hz. Via equation (2.6), the fundamental period at each element of the array was 


1 millisecond. The sampling parameter, S, was set equal to 2, 4, and 6, which 
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corresponded to a sampling frequency equal to the Nyquist Rate, twice the Nyquist 
Rate, and three times the Nyquist Rate, respectively. The number of time samples 


taken per element of the array is given by 
i, = Gieesl feu) 


where L is the total number of samples taken per element, A is the total number 
of harmonics present in the signal, and S is the sampling parameter. 

Baseline, or “no noise” test case results, were the first generated to ensure that 
the algorithm was working properly in a noise-free environment. Figure 3.1 depicts 
the noise free, time domain, received signal for case 2. For each baseline test case, 
and for a given value of the sampling parameter, bearing and depression angle 
estimation errors (measured in degrees) were obtained by running the computer 
simulation once and allowing the modified complex LMS algorithm 100 iterations. 
Case 1 baseline test case results showed zero degree estimation errors whereas case 
2 and case 3 baseline results showed estimation errors of less than 0.1 degrees. Case 
4 baseline results will be discussed later in this chapter. 

Following compilation of the baseline results, identical test cases were run 
using additive, wide-sense stationary, zero mean, white, guassian noise samples to 
corrupt the time samples of the received signal. Figure 3.2 depicts the time domain, 
received signal for case 2 fora SNR =O dB. For each test case, and for a given 
value of the sampling parameter and input SNR at a single element of the array, 
average bearing and depression angle estimation errors were obtained by running 
the computer simulation 100 times per SNA value and allowing the complex LMS 
algorithm 100 iterations per run. In all test cases, allowing the adaptive algorithm 


more than 100 iterations did not decrease the estimation errors. 
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Figure 3.1. Real Received Signal at Element (1,1), for Case 2 with No 
Noise Present. 
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Figure 3.2. Real Received Signal at Element (1,1) for Case 2 for SNR =0 
dqB. 
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Case 1 placed a single target at broadside relative to the planar array ( i.e., 
6, = 0°). This was considered the simplest case for the algorithm since it is 
at broadside that the far-field beam pattern has its narrowest beamwidth and, 
as a result, the algorithm should provide good angular resolution [Ref. 8: pp. 
62-66]. The general plane-wave field radiated by the target consisted of K = 6 
harmonics. Therefore, with AK = 6 and S = 2,4, and 6, only LZ = 13, 25, and 37 
time samples, respectively, were taken per element in the array (equation (3.1)). 
Figure 3.3 presents the average estimation error of the depression angle versus the 
sampling parameter S for the 0 dB SNR case. It can be seen that as the value of 
S increases for a given harmonic, the magnitude of the estimation error decreases. 
If S is held constant, it can be seen that as the harmonic number (q) increases, the 
magnitude of the estimation error decreases. These two trends were expected. As 
the value of S increases, more time samples are being processed and, as a result, 
the noise component of the received signal should tend to average out to be Zero. 
Since more data is available for evaluation, a better estimate of the deterministic 
signal results. An increase in the harmonic number, q, represents an increase in the 
frequency and, as a result, the beamwidth of the far-field beam pattern decreases 
(Ref. 8: pp. 58-62]. The decreased beamwidth increases the angular resolution of 
the algorithm and, as a result, decreases the estimation error. 

Figure 3.4 presents the results for case 1 at 9 dB SNR. The trends which 
were present in the 0 dB case are again apparent here. It also should be noted 
that as the SNR value increases (i.e., from 0 dB to 9 dB), the magnitude of the 
estimation error decreases. 

Results for bearing angle estimation error in the broadside case are irrelevant 


since the target is directly above the array, and the bearing angle has no meaning. 
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Figure 3.3. Average Depression Angle Estimation Error vs Sampling 
Parameter S, Case 1: SNR =0dB. I= 100. 
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Case 2 placed three broadband targets at random locations. The general plane- 
wave field radiated by each of the targets contained two unique harmonics. Target 
1 was located at (0, = 49°, = 38°) and radiated harmonics 1 and 6. Target 
2 was located at (9, = 5,y%, = 137°) and radiated harmonics 2 and 5. Target 
3 was located at (0, = 77°, = 307°) and radiated harmonics 3 and 4. Since 
there are three incident plane-wave fields each containing two unique harmonics, 
the output electrical signal from each element in the array exhibits a total of K = 6 
harmonics. As before, equation (3.1) dictates that only L = 13,25, and 37 time 
samples are to be taken per element of the array when S = 2,4, and 6, respectively. 
Figure 3.5 presents case 2 results for the depression angle estimation error with a 
0dB SNR. In general, as S increases, the estimation error decreases for a given 
harmonic. When comparing the harmonics associated with a particular target, 
the highest harmonic usually is associated with the lessor of the estimation errors. 
These general results can be explained using the same arguments presented for case 
1. 

Figure 3.6 presents the case 2 bearing angle estimation errors forO dB SNR. 
Figures 3.7 and 3.8 present the case 2 estimation errors for the depression and 
bearing angle, respectively, for 9 dB SNR. The same general trends already in- 
troduced are apparent in Figures 3.5 through 3.8. As in case 1, as the SNR value 
increased, the magnitude of the estimation errors decreased. 

Case 3 placed a single broadband target at endfire relative to the planar array 
(0, = 90°, = 90°). This case was considered the most difficult for the algorithm 
since it is at endfire that the far-field beam pattern beamwidth is its broadest [Ref. 
8: pp. 62-66]. The general plane-wave field radiated by the target consisted of g = 5 
harmonics. Therefore, with S = 2,4, and 6, only LD = 11,21, and 31 time samples 


were taken, respectively, at each element of the array (equation (3.1)). Figures 
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28 


45.0 


: 

LEGEND 
j a gr 
2g O = g=2 
ww =) A= ae 
> an Gime 
5 : 2 es 
L- © a= 9=G 
eA 
O 
oO 
Oo 
iJ 
= 
@ 
—~ © 
mw 
q- 
= 7 
a A 
ee 
(ti) 
(5 2 
co; 
& 
> os 
c a 





2.0 4.0 2) 510 
ST aeNG Senin cle VEioee NEOUS Tiss ) 


Figure 3.6. Average Bearing Angle Estimation Error vs Sampling Pa- 
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Figure 3.8. Average Bearing Angle Estimation Error vs Sampling Pa- 
rameter S, Case 2: SNR=9dB I=100. 
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3.9 and 3.10 illustrate the average depression and bearing angle estimation errors, 
respectively, for 0 dB SNR. Figures 3.11 and 3.12 illustrate the same information 
for9 dB SNR. The estimation errors for the highest harmonic (q = 5) are not 
shown for any of the case 3 results since the magnitude of these errors was between 
70 and 90 degrees. This poor performance for the highest harmonic at endfire can 
be explained theoretically (see Chapter 2 - phase unwrapping). 

Case 4 corresponded to three broadband targets being present. The general 
plane-wave field radiated by each target was composed of two harmonics. In this 
case, two of the targets shared a common spectral line, namely, the fundamental 
frequency (¢ = 1). Target 1 was located at (6, = 45°,y~, = 0°) and radiated 
harmonics 1 and 2. Target 2 was located at (0, = 45°,~, = 180°) and radiated 
harmonics 1 and 5. Target 3 was located at (@, = 33°,~, = 47°) and radiated 
harmonics 3 and 4. A total of A = 5 harmonics were areecnt in the output 
electrical signals, and, with S = 2,4, and 6; only 2 = 11,21, and 31 time samples, 
respectively, were taken at each element of the array (equation (3.1)). For this 
case, only no noise results were compiled, with the computer simulation running 
once and the complex LMS adaptive algorithm allowed 100 iterations. 

For harmonics g = 2,3,4, and 5 contained in the output electrical signal, 
the algorithm correctly identified the location of the targets with zero degrees 
estimation error. However, for the shared harmonic (q = 1), the algorithm located 
a false fourth target, which was located exactly between the actual locations of the 
targets whose radiated general plane-wave fields contained the shared harmonic, 
namely, targets 1 and 2. 

Tabular numerical data for test cases 1 through 4 is contained in The Ap- 


pendix. 
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Figure 3.9. Average Depression Angle Estimation Error vs Sampling 
Parameter S, Case 3: SNR=0dB I[=100. 
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Figure 3.10. Average Bearing Angle Estimation Error vs Sampling Pa- 
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Figure 3.11. Average Depression Angle Estimation Error vs Sampling 
Parameter S, Case 3: SNR=9dB I= 100. 
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Figure 3.12. Average Bearing Angle Estimation Error vs Sampling Pa- 
rameter S, Case 3: SNR=9dB [= 100. 
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test case 2 because it was felt that this case presented the most realistic real world 
situation. First, a minimum number of iterations for the algorithm to produce the 
same results as presented in this chapter (when 100 iterations of the algorithm was 
allowed) was found. The minimum number was found to be 10 iterations. Using 
this information, the step size parameter, yu, (which had been a constant) was 
varied using a step function. For the first 10 iterations of the algorithm, p was held 
constant. This allowed the algorithm to reach its minimum value with a constant 
yw. The value for » was then stepped down to one-tenth its original value for an 
additional 90 iterations. No decrease in estimation error resulted. The parameter 
was then stepped down to one one-hundredth, and then one one-thousandth, of 
its original value, with no change in the estimation error resulting. Additional 
attempts to model the step size parameter as a decreasing exponential yielded no 


change in the estimation error. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


The goal of this thesis was to evaluate a frequency domain adaptive beamform- 
ing algorithm which was developed by Ziomek and Chan [Ref. 5]. The multiple 
broadband target localization capability and the full angular coverage capability 
of the algorithm were of particular interest. 

Chapter 3 presented the results of four test cases which were designed to 
test the algorithm for several different capabilities. Several conclusions concerning 
the performance of the frequency domain adaptive beamforming algorithm can be 
drawn from these results. Among these conclusions is the fact that the algorithm 
performs well in all test cases, yielding what would be considered acceptable esti- 
mation errors for the target localization problem. It should also be noted here that 
the number of time samples taken per element of the array never exceeded 37 time 
samples, an amount which represents a small number of data points. 

The algorithm does exhibit multiple broadband target localization capability. 
The test case results show that multiple targets can be localized if their radiated 
acoustic plane-wave fields contain at least one unique spectral line. The separation 
distance of the spectral lines can be controlled by selecting an appropriate value for 
the FFT bin spacing. The algorithm has also demonstrated its ability to localize 
targets both at endfire and broadside positions relative to the planar array, thus 
demonstrating the algorithm’s full angular coverage capability. 

Several general trends were apparent in the results. The first of these trends is 
that as the sampling parameter S increased, the magnitude of the estimation error 
decreased. This was due to a corresponding increase in the amount of time samples 


being processed. The second trend observed was that as the harmonic number q 
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increased, the magnitude of the estimation error, in general, decreased. This trend 
can be explained by the fact that as the harmonic number q increases (equivalent to 
an increase in the frequency), the far-field beam pattern beamwidth decreases which 
results in better angular resolution. The third trend observed was that asthe SNR 
increased (the output electrical signal becoming less noisy), the magnitude of the 
estimation error decreased. It was also found that the estimation errors presented 
in this thesis are obtainable after ten iterations of the adaptive algorithm, and an 
increase in the number of iterations did not decrease the estimation error. 

In the course of this investigation, several possible areas for future research 


presented themselves: 


e further study of the effects of varying the step-size parameter, pz, in an attempt 
to decrease the magnitude of the estimation error, 


e investigation of a noise reduction system prior to processing the output elec- 
trical signals, and 


e application of other spectral analysis techniques (i.e., autoregressive, maxi- 
mum entropy, maximum likelihood, etc.) to produce frequency spectra. 


39 


APPENDIX 


The Appendix presents the tabular numerical data from which the graphical 
results of Chapter 3 were produced. The numerical data for the graphical results 
of Figures 3.3 through 3.12 is presented in Tables A.1 through A.10, respectively. 
Tables A.11 and A.12 presents the numerical data for Case 4. There is no corre- 


sponding graphical result for Case 4. 
TABLE A.1. NUMERICAL DATA CORRESPONDING TO FIGURE 3.3. 


Case 1: One target located at broadside relative to planar 
array; Six harmonics present in output electrical signal; 
sampling parameter S = 2, 4, and 6; 100 iterations; SNR = 0dB 


HARMONIC LOCATION DEPRESSION ANGLE 6, 
aR ESTIMATION ERROR (DEG) 
(DEG) (DEG) S=2 S=4 S=6 


270 -7.388 -4732 -4.124 
2/0 23.207 Meee 0 Selo 
Zo “2.542 -1654 -1.374 
2/70 =2 50S 225) ano oS 
270) -1.312 -0.962 -0.784 
Ze =| 086350592 —O1660 
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TABLE A.2. NUMERICAL DATA CORRESPONDING TO FIGURE 3.4. 


Case |: One target located at broadside relative to planar 
array, SIX harmonics present in output electrical signal; 
Sampling parameter S = 2, 4, and 6; 100 iterations; SNR = 9dB 


HARMONIC LOCATION Der eosl ON) FINGEE Ye: 
Oo, v, ESTIMATION ERROR (DEG) 
(DEG) (DEG) SZ 9=4 SS 


ZO eee 2 a1. /90,>1\2456 
ZO -1.191 -Q0.94Q0 -0.686 
2/0 —Oy69° ~=0562 —0.449 
210 -0.594 -0.446 -0.365 
2/0 =O BeO.526 --0 207 
270 Or OC) ie, 516-0724 





TABLE A.3. NUMERICAL DATA CORRESPONDING TO FIGURE 3.5. 


S2ce 2)  Ultiple targets located at Fandom positions; 
Six nechmlem! cS PrESenl in GUEDUL Electrical signal; 
oampling parameter 5 = 2, 4, and 6; 100 iterations; SNR = OdB 


HARMONIC LOCATION DEPRESSION ANGLE 6, 
aA 4v, ESTIMATION ERROR (DEG) 
(DEG) (DEG) S=2 S=4- S=6 


49 38 a en) (OS --0.25s 
» i S7 OOOO 450  —-O,522 
vy SOY Tigo o.00)] . 1.075 
Le BOY ino lemme |,OG5 G52) 
5 bow -0.599 -0.164 -0.010 
49 ENS Zee DL OOO 
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TABLE A.4. NUMERICAL DATA CORRESPONDING TO FIGURE 3.6. 


Case 2: Multiple targets located at random positions: 
SIx harmonics present in output electrical Signal; 
sampling parameter S = 2, 4, and 6; 100 iterations: SNR = OdB 


HARMONIC EOCHIIGN BEARING ANGLE &, 
eh b, ESTIMATION ERROR (DEG) 
(DEG) (DEG) Sa 9=4 Sas 


49 38 nOWZO | O20 omc. 00) 
» 137 OR 7 San 2 Smal O77 
V7 5O7 +43.268 +27662 +14599 
U7 S017 $2020 5.769 17100 
=) 17 HOAGS el Wee Or Or 
49 38 =O1420° 38220) +0207 





TABLE A.5. NUMERICAL DATA CORRESPONDING TO FIGURE 3.7. 


Case 2: Multiple targets located at random positions; 
Six harmonics present in output electrical signal; 
sampling parameter S = 2, 4, and 6; 100 iterations; SNR = 9dB 


HARMONIC LOCATION DEPRESS! ON Ghee: 
Ss b, ESTIMATION ERROR (DEG) 
(DEG) (DEG) iV 9=4 S16 


AQ 38 -0.338 -0.113 -0.080 
5 Ss -0.053 -0.073 -0.050 
Li SOF -0.425 +0280 -0.192 
77? 3507 +1.383 +0.372 +1.350 
5 [Se +0.017 -0.013 -0.003 
AQ 38 -0.447 +0.059 +0.022 
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TABLE A.6. NUMERICAL DATA CORRESPONDING TO FIGURE 3.8. 


Case 2 Multiple targets located at random positions: 
Six harmonics present in output electrical signal; 
Sampling parameter S = 2, 4, and 6; 100 iterations: SNR = 9dB 


HARMONIC LOCATION BEARING ANGLE y, 


© 


.°] 


Yo 


(DEG) (DEG) 


49 
5 
a) 
Ua 
S 
49 


38 
| 97 
BO 
DO 
7 
38 


ESTIMATION ERROR (DEG) 


SS 2 


WZ 
2 || Aexe: 
SOO 46 
=e 2 | 
+0.648 
AON OSS: 


S=4 


+0.048 
mOZ69 
CORY 
50.075 
-0.148 
FOZ 


S18 


-Q0.145 
Oe) 
BOLO SO 
alle oo 
OS eke 
OL OS 





TABLE A.7. NUMERICAL DATA CORRESPONDING TO FIGURE 3.9. 


aoe 5. One target located at endiire nelative to planar array; 
five harmonics present in output electrical signal; 
Pameuna parameter >= 2,4, and 6; lOO iterations; SNR = OdB 


LOCATION DeeReostON ANGLE 9, 
sl. b, ESTIMATION ERROR (DEG) 
DEG) (Dac) SZ S=4 50 


HARMONIC 


Soli ++ 72559 
+15.058 +8.046 
46.119 +4723 
fy o4 “F2,.926 
Oo, 709470 85 | 


eH Oke: 
tie 7S 


90 90 
90 90 
90 +7594 
90 +15.749 
90 Ce) Z 
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TABLE A.8. NUMERICAL DATA CORRESPONDING TO FIGURE 3.10. 


Case 3: One target located at endfire relative to planar array; 
rive harmonics present in output electrical signal; 
sampling parameter S = 2, 4, and 6; 100 iterations; SNR = OdB 


HARMONIC LOCATION BEARING ANGLE &, 


Q v, ESTIMATION ERROR (DEG) 
(DEG) (DEG) S=2 S=4 S=6 


90 90 
90 90 
90 90 
90 90 
90 90 


-0.547 
= EZ8Z 
SOmOS 
mee) asl 
-69.48 | 


= 0) 5) 7 Sg 0) 21), 25 
Ono. 7 Zan 
-0.038 -0.094 
“Oia -0.04i 
SODCIS aso 587 





TABLE A.9. NUMERICAL DATA CORRESPONDING TO FIGURE 3.11 


Case 3: One target located at endfire relative to planar array; 
five harmonics present in output electrical signal; 
sampling parameter 5S = 2, 4, and 6; 100 iterations; SNR = 9dB 
HARMONIC LOCATION DEPRESSION ANGLE Ge: 
ee bs, ESTIMATION ERROR (DEG) 
(DEG BEas) 9=2 S=4 S=6 


90 90 
90 
90 
90 
90 


23) 2)(C!] 
ae OS) 20) 
+3.744 
+2.483 


+68.938 


+4994 
se Sel AS 
tZO2 7 
eZ 5972 


+66.092 


+4677 
+2.932 
5647/0) | 
‘Zo 2 
*7Oiea0 





TABLE A.10. NUMERICAL DATA CORRESPONDING TO FIGURE 3.12. 


Case 5: One target located at endfire relative to planar array; 
rive harmonics present in output electrical signal; 
sampling parameter S = 2, 4, and 6; 100 iterations; SNR = 9dB 


HARMONIC LOCATION 


Go 8, 
(DEG) (DEG) 


90 90 
90 90 
90 90 
90 90 
90 90 





BEARING ANGLE ¥&, 
ESTIMATION ERROR (DEG) 
S=2 S=4 S=6 


0:0. O9/ =-0 0S 
“JIOSO) (0) | Beis 510) (Osis 
-0.090 -0.045 +0019 
HOWSe + O009° —0,021 


= (OG Oiz e890 250 595,056 


TABLE A.11. NUMERICAL DATA FOR DEPRESSION ANGLE ESTIMATION 
ERRORS FOR CASE 4. 


Case 4: Multiple targets; targets one and two share a common 
Seecural lime; five harmonics present in Output electrical signal; 
Samelingparameter > = 2, 4 and 6; 100 iterations; no noise 


HARMONIC LOCATION 


Oo =, 
(DEG) (DEG) 


45 QO 
45 0 
45 
45 
5S 
SS) 





DER RGoslON FINGEE “GF 
ESTIMATION ERROR (DEG) 
SEZ 5-4 SSS 


ml Oooo) 9605 
0.0 0.0 0.0 
TICs ee OSS 5+ |9 665 
0.0 0.0 0.0 
0.0 0.0 0.0 
0.0 0.0 0.0 


TABLE A.12. NUMERICAL DATA FOR BEARING ANGLE ESTIMATION 
ERRORS FOR CASE 4. 


Case 4 Multiple targets; targets one and two share a common 
spectral line; five harmonics present in output electrical signal; | 
Sampling parameter 5 = 2, 4, and 6; 100 iterations; no noise 


HARMONIC LOCATION BEARING ANGLE ¥, 
6 ov, ESTIMATION ERROR (DEG) 
(DEG) (DEG) $=2 S=4  S=6 


45 -| 60100 -ls0roo: — leo 
45 0.0 0.0 
45 0.0 0.0 
45 0.0 0.0 
5S 0.0 0.0 
SS 0.0 0.0 
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